Bohr Hamiltonian with deformation-dependent mass term for the Davidson potential 
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Analytical expressions for spectra and wave functions are derived for a Bohr Hamiltonian, de- 
scribing the collective motion of deformed nuclei, in which the mass is allowed to depend on the 
nuclear deformation. Solutions are obtained for separable potentials consisting of a Davidson poten- 
tial in the ft variable, in the cases of 7-unstable nuclei, axially symmetric prolate deformed nuclei, 
and triaxial nuclei, implementing the usual approximations in each case. The solution, called the 
Deformation Dependent Mass (DDM) Davidson model, is achieved by using techniques of supersym- 
metric quantum mechanics (SUSYQM), involving a deformed shape invariance condition. Spectra 
and B(E2) transition rates are compared to experimental data. The dependence of the mass on 
the deformation, dictated by SUSYQM for the potential used, reduces the rate of increase of the 
moment of inertia with deformation, removing a main drawback of the model. 



I. INTRODUCTION 

The Bohr Hamiltonian [T] and its extensions, the geo- 
metrical collective model [2j [3] , have provided for several 
decades a sound framework for understanding the collec- 
tive behaviour of atomic nuclei. It has been customary 
to consider in the Bohr Hamiltonian the mass to be a 
constant. However, evidence has been accumulating that 
this approximation might be inadequate. In particular: 

1) The moments of inertia are predicted to increase 
proportionally to /3 2 , where j3 is the collective variable 
corresponding to nuclear deformation, while the experi- 
mentally determined (from the spectra) moment of iner- 
tia shows a much more moderate increase as a function of 
the experimentally determined (from the B(E2) transi- 
tion rates) deformation, especially for well deformed nu- 
clei [1]. This discrepancy has led to arguments that the 
use of the Bohr Hamiltonian is justified for vibrational 
and transitional nuclei, but its applicability to deformed 
nuclei needs further clarification. 

2) Detailed comparisons to experimental data have re- 
cently pointed out [SJ [5] that the mass tensor of the col- 
lective Hamiltonian cannot be considered as a constant 
and should be taken as a function of the collective coordi- 
nates, with quadrupole and hexadecapole terms present 
in addition to the monopole one. 

3) In the framework of the Interacting Boson Model 
(IBM) 7] , which offers an algebraic description of atomic 
nuclei complementary to that of the Bohr Hamiltonian, 
it is known that in its geometrical limit [7. , obtained 
through the use of coherent states [7] , terms of the form 
/3 2 7r 2 and/or more complicated terms appear [8], in addi- 
tion to the usual term of the kinetic energy, 7r 2 . Thus it 
might be appropriate to search for a modified form of the 
Bohr Hamiltonian, in which the kinetic energy term will 
be modified by terms containing /3 2 and/or more compli- 
cated terms. 

Based on this evidence, a Bohr Hamiltonian with a 



mass depending on the collective variable (3 can be con- 
sidered. Position-dependent effective masses have been 
studied recently in a general framework [9] , while several 
Hamiltonians known to be soluble through techniques 
of supersymmetric quantum mechanics (SUSYQM) [101 
ITT] , have been appropriately generalized [12] to include 
position-dependent effective masses, the 3-dimensional 
harmonic oscillator being among them |12) . 

In the present work we are going to show that a Bohr 
Hamiltonian with a Davidson potential |13j in j3 (a har- 
monic oscillator potential with a term proportional to 
1//3 2 added to it) can be generalized in order to include 
a mass depending on /3, B = B /(l + a/3 2 ) 2 , where B 
and a are constants. We shall call this approach the 
Deformation Dependent Mass (DDM) Davidson model. 
Three cases of potentials, for which exact separation of 
variables can be achieved, will be considered: 

a) Potentials independent [H] of the collective variable 
7 (an angle measuring departure from axial symmetry), 
called 7-unstable potentials, appropriate for describing 
vibrational and near-vibrational nuclei. 

b) Potentials of the form [HHT8] v(@,j) = u(j3) + 
w(j)/j3 2 , with u{(3) being the Davidson potential [15] , 
and with w("f) having a deep minimum at 7 = 0, corre- 
sponding to axially symmetric prolate deformed nuclei. 

c) Potentials of the form v(f3, 7) = u{(3) + 10(7)/ '/3 2 , with 
u((3) being the Davidson potential [13], and with 10(7) 
having a deep minimum at 7 = 7r/6, corresponding to 
triaxial nuclei [TH1 [20] ■ 

Analytical results for spectra and B(E2) transition 
rates will be provided for all three cases, implement- 
ing the usual approximations in each limit |21H23] , while 
comparison to experimental results will be undertaken in 
the first two, for which able bulk of experimental data 
exists. A special solution regarding 7-unstable nuclei has 
been given earlier in Ref. [2~i] . 

The analytical spectra and wave functions of the Bohr 
Hamiltonians considered are obtained by using tech- 
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niques of supersymmetric quantum mechanics jTDl [TTj . 
equivalent [EQ to the factorization method of Infeld 
and Hull [25] ■ The integrability of the Hamiltonian is 
achieved by imposing a deformed shape invariance con- 
dition [T5]. These tools are described in more detail in 
Section VI. 

It should be noticed that the concept of a non-constant 
mass in the framework of the Bohr Hamiltonian has been 
used long ago in numerical solutions of a generalized Bohr 
Hamiltonian [26], as well as in relevant mean field calcu- 
lations [27 . The main difference of the present work from 
these earlier approaches is that analytical solutions are 
obtained here. In addition, in the present case the num- 
ber of free parameters remains small (two or three) , while 
the functional dependence of the mass on the deforma- 
tion for the potential used is dictated by SUSYQM. The 
relation of the present work to these earlier approaches 
will be discussed in Section XII. 

The structure of the present work is as follows. In 
Section II the formalism of position-dependent effective 
masses, which we use in order to allow the mass to depend 
on the deformation /?, is briefly reviewed, and applied to 
the Bohr Hamiltonian in Section III. The three exactly 
separable cases described above are considered in Sec- 
tion IV, in which the common overall form of the radial 
equation in all three cases is pointed out, while in Sec- 
tion V we focus on the use of the Davidson potential in 
the radial equation. The solvability of the Hamiltonian is 
achieved in Section VI by imposing a deformed shape in- 
variance condition, leading to the energy spectrum given 
in Section VII and the wave functions given in Section 
VIII. Normalization coefficients are given in Section IX, 
while a detail on their numerical calculation is included 
as Appendix 1. B(E2) transition probabilities are con- 
sidered in Section X, while in Section XI comparisons 
of spectra and B(E2)s to experimental data are carried 
out. Finally, connections to earlier work are discussed in 
Section XII, while Section XIII contains discussion of the 
present results and plans for further work. 



Hamiltonian 



H 2 

H = - — [m 5 (x)Vm K (x)Vm A (x) 
(x)Vin"'(x)Vm f (x)] + V(x), 



(1) 



where V is the relevant potential and the parameters 5' , 
k', A' are constrained by the condition 6' + k' + A' = — 1. 
Assuming a position dependent mass of the form 



m(x) = mo-M(x), M(x) 



1 



(/(x)) 5 



/(x) = l + . 9 (x), 
(2) 

where m is a constant mass and M(x) is a dimensionless 
position-dependent mass, the Hamiltonian becomes 



H = - — [/ 5 (x)Vr(x)V/ A (x) 
4m 

+/ A (x)V/ K (x)V/ 5 (x)] + V(x), 



(3) 



with 5 + k + A = 2. It is known [9] that this Hamiltonian 
can be put into the form 



H = - 



2mo 



77(xTv/(x)vy/(x) + y e// (x), (4) 



with 



V e// (x) = V(x) + 



s>2 r 



2r?io 



-(1-<5-A)/(x)V 2 /(x) 



\-6 



- - A ] (V/(x)) 2 



(5) 



where 5 and A are free parameters. 

In the final part of the paper, in which comparison to 
experiment will be carried out by fitting the theoretical 
predictions to the experimental data, it will be seen that 
the predictions for the theoretical spectra turn out to be 
independent of the choice made for 8 and A. 



II. FORMALISM OF POSITION-DEPENDENT 
EFFECTIVE MASSES 



For reasons of completeness, we briefly review the ba- 
sics of the formalism needed in handling effective masses 
depending on the position. The main problem encoun- 
tered is the generalization of the kinetic energy term. We 
show how this can be solved in an unambiguous way. 

When the mass m(x) is position dependent [9], it does 
not commute with the momentum p = —iKV. There- 
fore, there are many ways to generalize the usual form 
of the kinetic energy, p 2 /(2mo), where mo is a constant 
mass, in order to obtain a Hcrmitian operator. In or- 
der to avoid any specific choices, one can use the general 
two-parameter form proposed by von Roos |28| . with a 



III. BOHR HAMILTONIAN WITH 
DEFORMATION-DEPENDENT EFFECTIVE 
MASS 



A. Deformation-dependent effective mass 
formalism 



The original Bohr Hamiltonian pQ is 



H B = - 



2B 



JL_J^ 4 d_ _JL d_ . 

WdP + /3 2 sin37C^ Sln 7 ^y 



4/3 



- y 



Ql 



fe=l,2.3 



sin 2 (7 — §7rfc) 



v{p,-r), (6) 



where j3 and 7 are the usual collective coordinates (/3 be- 
ing a deformation coordinate measuring departure from 
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spherical shape, and 7 being an angle measuring depar- 
ture from axial symmetry), while Qk (k = 1, 2, 3) are the 
components of angular momentum in the intrinsic frame, 
and B is the mass parameter, which is usually considered 
constant. 

We wish to construct a Bohr equation with a mass de- 
pending on the deformation coordinate /?, in accordance 
with the formalism described above, 



B{P) = 



Bp 

(/(/3)) 2 



(7) 



where Bq is a constant. We then need the usual Pauli- 
Podolsky prescription [53] 



(V3)' = g*> 



<9$ 



V 2 $ = -Ldiy/g^d^, (8) 
\J9 



in order to construct a Schrodinger equation correspond- 
ing to the Hamiltonian of Eq. Q in a 5-dimensional 
space equipped with the Bohr- Wheeler coordinates /3, 7. 
Since the deformation function / depends only on the ra- 
dial coordinate /3, only the (3 part of the resulting equa- 
tion will be affected, the final result reading 



P 
+ - — 

8^ 2 



d f 2 d d 

2 /3 4 8(3 ^ f 0/3 ^ ~ 2/3 2 sin 3 7 dj ^ ^ dl 



2 

E 



fe=l,2,3 



sin 2 (7 — §7rfc) 



"eff 



(9) 



where reduced energies e = BqE/H 2 and reduced poten- 
tials v = BoV/h 2 have been used, with 



where 

ds 2 = gijdqidqj, 
the symmetric matrix 3^ having the form 



(9ij) 



/5n 9i2 513 \ 

921 922 

531 g 33 

344 

V 555 / 



with 130] 





_ Jl 


511 


B 






512 


-s< 




_ Ja 


513 


" B 




_ Ji 


522 


B 




_ Js 


533 


B 


544 


= 1, 


555 


= /3 2 , 



Ji 



= ^ sin^ 6 cos^ ip + — sin 2 sin^ ^1 



cos 6, 



J2 



sin -0 + -ft cos V 



where the moments of inertia are 



(13) 



(14) 



Jz _ j 



cos 9, 



(15) 



"e// = VW, l) + -Al-5-\)fV 2 f 



A (V/) 2 



(10) 



J k = AB(3 2 sin 2 7 - k 



2tt 



The determinant of the matrix is 



(16) 



B. Connection to curved space 

In Ref. [9] it has been proved that the position- 
dependent effective mass formalism can be equivalently 
expressed in a curved space. We shall prove here that 
this connection is possible also in the case of the Bohr 
Hamiltonian, paving the way for connecting in Section 
XII the present results to earlier related work. 

Ordering the coordinates as 

qi = q_2 = 0, 33 =ip, 54 = P, <75 = 1- (11) 

the kinetic energy in the standard Bohr Hamiltonian [T] 
can be represented as 



T = 



B I 'ds 
dt 



(12) 



g = Jl ^ 3 Ja p 2 sin 2 8 = 4/3 8 sin 2 3 7 sin 2 6. (17) 



The relevant volume element is then 



dV = 2/3 4 sin 37 sin Qd^dQd^dPd-y. 



The inverse matrix is found to be 



(91 



(9-11 1 


512 1 


513 1 


921 


5 22 J 


523 1 


531 1 


532 1 


533 1 












V 



\ 





744 1 

555V 



(18) 



(19) 
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with 



9U 1 



9x2 



9x3 



922 = 



923 



B 



cos 2 ip sin ip 



sin 2 B V Ji J2 

( 1 X \ sin ip cos ip 

B 

-B 



Jx J 2 J sin6 
cos 2 ip t sin 2 ip \ cot 8 
Jx 



J 2 J sin9' 



sin ip cos 2 ip 



Ji 

1 1 

J~X ~ J~2 



J2 

cot @ sin ip cos ip, 



(20) 



533 



9il 

955 



_ / cos 2 ip sin 2 i/> \ 2 „ B 

1, 

1 

is 5- 



Using these matrix elements and the value of the deter- 
minant from Eq. (17) in Eq. (fsl) we obtain 



T 



2B 

1 d 



2B 

a 



/3 4 <9/T d(3 



Ql 



/3 2 sin 3 7 d 7 7 S 7 4/3 2 , ^ Q sin 2 (7 - f Trfc) 

fc— 1,2,0 ^ " ' 



(21) 



where are the components of the angular momentum 
in the intrinsic frame 



Qx i 

Qy ^ 

Qz = -i 



cos ip d , d „ .9 

"^~^^ + sm ^:SH +cotecos?A— , , 
sin B a$ c?B aV> / 

sin 7/1 9 ,9 . 9 

— ;7~ht: + cos ^ a?^ - cote sin 
sm B a$ aB dip 



dip ' 



(22) 



The connection between the position-dependent effec- 
tive mass and curved spaces has been considered in Ref. 
[5]. According to the findings of Ref. [5], one expects 
in the present case all elements of the matrix ( 14 ) to be 
divided by / 2 



q>.. = Hi 

As a result, the determinant of the matrix will be 

' = JL 
9 f xo> 

and the volume element will be 



(23) 



(24) 



(25) 



The elements of the inverse matrix will be 
1 t2„-x 



9 



1 .1 



f A 9 



(26) 



According to Ref. 9 , in order to obtain the 
Schrodinger equation in the form of Eq. ([9]), one has 
to start with the equation 



H g * = 



where 



^ = 



1 1 

2V.9' 



* = / 5/2 *, 



(27) 
(28) 



while reduced energies and reduced potentials are used, 
as in Eq. Q. The exponent in the last equation is related 
to the dimensionality of the space. 

Substituting the g' matrix elements and determinant in 



Eq. (27), and performing the relevant calculation (which 
closer 



y resembles the pure Bohr case, except for the 44- 
term), we see that Eqs. (27) and ^ do coincide with 



/ = 



df 



f = 



d 2 J 



u g =u eff + ff- 2(/)»+4= . 

This result has several important consequences. 

1) It becomes clear that solving the Schrodinger equa- 
tion ([9| with deformation dependent mass is equivalent 
to solving a modified Bohr equation (27) with different 
metric matrix g' and another effective potential, u g . Be- 
tween the two equivalent schemes, one chooses to solve 
Eq. ([9]) instead of Eq. (27), just because the former 
can be solved analytically through the use of SUSYQM 
techniques. 

2) The wave functions \& = J 5 / 2 "]/ are accompanied by 
the volume element dV — dV/ f 5 . As a result 



J %*VdV' = J (/ 5 / 2 **)(/ 5 / 2 * 



dV 

IF 



= / ^*^dV 7 



(30) 

i.e., the wave functions \& of the deformation dependent 
mass problem correspond to the usual Bohr volume ele- 
ment dV. 

3) The simple relation between ^ and "J/ also shows 
that the wave functions 'J satisfy the well-known 24 sym- 
metries of Bohr wave functions p] , which the wave func- 
tions "J satisfy by construction. If these symmetries were 
not satisfied, the solutions could not have been used for 
the description of nuclei. 

Further consequences, regarding the connection of the 
present approach to earlier work, will be discussed in 
Section XII. 



IV. EXACTLY SEPARABLE SPECIAL FORMS 
OF THE BOHR HAMILTONIAN 

The solution of the above Bohr-like equation can be 
reached for certain classes of potentials using techniques 
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developed in the context of SUSYQM [T0HI3]. At this 
point exact separation of variables can be achieved in 
three cases, described in the following three subsections. 



A. 7-unstable nuclei 

In order to achieve separation of variables we assume 
that the potential 7) depends only on the variable 
j3, i.e. v(j3) = u{j3) [14] . Potentials of this kind are 
called 7-unstable potentials, since they are appropriate 
for the description of nuclei which can depart from axial 
symmetry without any energy cost. 

One then seeks wave functions of the form [T3J [3T] 



*C8,7^ i )=CC8)*(7 ) fli) ) 



(31) 



with 111(7) having a deep minimum at 7 = 0. Then the 
angular momentum term can be written as |21j 



;(Ql + Q 



E 

fc=l,2,3 

Ql) + Ql 



Qt 



■ 2 / 

sin (7 - 



sin 2 7 



§**) 



One then seeks wave functions of the form [! 

W7^) = <^(/3,7)2?M./f(^), 



(35) 



(36) 



where T>(9i) denote Wigner functions of the Euler angles, 
L is the angular momentum quantum number, while M 
and K are the quantum numbers of the projections of 
angular momentum on the laboratory-fixed z-axis and 
the body-fixed z'-axis respectively. Then separation of 
variables leads to 



where 0t (i — 1, 2, 3) are the Euler angles. Separation of 
variables gives 



2 /3 4 df3 P J 8(3 



- S 



- A 



P 

+ - — 
2/3 2 

(V/) 2 



A 



6-X)fV 2 f 
(32) 



1 d a 

sm 37 07 07 



4 ^ sin 2 



$(7,1 



_Ql 

(7 - ^k) 

)=A$( 7 ,^) 



(33) 



Eq. ( 33 ) has been solved by Bes 



A = t(t + 3) rep- 
resents the eigenvalues of the second order Casimir oper- 
ator of SO (5), while r is the seniority quantum number, 
characterizing the irreducible representations of SO (5). 
The values of angular momentum L occurring for each 
r are provided by a well known algorithm and are listed 
in 0Q3]. Within the ground state band (gsb) one has 
L = 2t. The L = 2 member of the quasi-71 band is de- 
generate with the L = 4 member of the gsb, the L = 3, 
4 members of the quasi-71 band are degenerate to the 
L = 6 member of the gsb, the L = 5, 6 members of the 
quasi-71 band are degenerate to the L = 8 member of 
the gsb, and so on. 



B. Axially symmetric prolate deformed nuclei 

In order to achieve exact separation of variables, we 
assume a potential of the form [TH - fTS] 



(34) 



2 /3 4 8(3 

4G-« 



1 



dp 



f 



1 



2 



■-— -— sin 37— 
sin 37 07 07 



2/3 2 
(V/) 2 



i^ 2 



efiCS), (37) 



sin 2 7 



-2^(7)]7 ?K ( 7 )=Ar ?i f(7), (38) 



where 



A = A 



L(L + l) 



(39) 



and (frxiPii) = Cl(P)vk(i)- We remark that Eq. Q 
has the same form as Eq. |32l, obtained in the case of 



7-unstable nuclei, when A in the former is replaced by A 
in the latter. However, the results are different as far as 
the physics described is concerned. The angular momen- 
tum dependence, contained in A and A respectively, is 
different. Furthermore, the angular equation is different 
in each case, due to the different treatment of the 7 vari- 
able, the potential being confined to 7 ~ in the former 
case, while being independent of 7 in the latter. 

Eq. ( 38 ) has been solved for a harmonic oscillator po- 
tential 



™(7)^(3c)V, 



(40) 



in the case of 7 ss [TH1 [2T] , resulting in 
K 2 

A = e 7 -— , e 7 = (6c)(rc T + l), n 7 = 0, 1, 2, . . . , 

(41) 

where n 7 is the quantum number related to 7-oscillations. 
The allowed bands are characterized by 



0, K = 0: 



1. 



K = ±2; 



K = 0,±4; 



(42) 
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As a result 

- L(L + 1)-K 2 < L(L + l)-if 2 

A= — — he 7 = h(6c)(n 7 + l). 

(43) 



C. Triaxial nuclei with 7 = tt/6 

In this case we assume again a potential of the form of 
Eq. (34), but with w("f) having a deep minimum at 7 = 



7r/6. In this case K, the angular momentum projection 
on the body-fixed z'-axis, is not a good quantum number 
any more, but a, the angular momentum projection on 
the body-fixed rr'-axis, is a good quantum number, as 
found 22J in the study of the triaxial rotator []j|] |2"0] . 
Then the angular momentum term can be written as [221 



where n 7 is the quantum number related to 7-oscillations. 
As a result 



A 



AL(L + 1) - 3a 2 



'2c 



(51) 



We remark that Eqs. (37) and (461 have the same 



form, with A in the former replaced by A in the latter. 

In the literature on triaxial nuclei it is customary, in- 
stead of the projection a of the angular momentum on 
the a/-axis, to introduce the wobbling quantum number 
[21 [22] n w = L — a. Inserting a = L — n w in Eq. (51 ) 
one obtains 



s L{L + A) + Zn w {2L-n w ) ( 1 

A = + V2c I n 7 + - 



D. Common form of the radial equation 



(52) 



fc=l,2,3 ' 



( 7 - §7Tfc) 



4(Qi + Q2 + O3) — 3Qi (44) ^ 6 remar ^ ^hat (32), (37), and (46) have the 



One then seeks wave functions of the form [35] 

$(AT,fi)=^(AT)^, a («i), (45) 

where X>(#i) denote Wigner functions of the Euler angles, 
L is the angular momentum quantum number, while M 
and a are the quantum numbers of the projections of 
angular momentum on the laboratory-fixed z-axis and 
the body-fixed x'-axis respectively. Then separation of 
variables leads to 



2 p 4 dp' J dp 



f 



1 



/+^A+-(l-5-A)/V 2 f 



20 



\{\-> 



~-A)(V/V »( >'> 



a,a(/?) 

^L,a(/3), (46) 



- . — sin 37— + 2^(7) 
sm 37 07 07 



with ^08,7) = ZlAPMi), and 



77(7) = ^77(7), (47) 



7 4L(L + l)-3a 2 A , 
A = — h A'. 



(48) 



Eq. ( 47 ) has been solved for a harmonic oscillator poten- 
tial 

W ( 7 ) = |c( 7 -^) 2 , (49) 
in the case of 7 « 7r/6 [23] , resulting in 

(50) 



A' = e 7 = V2c ( n 7 , ^ . 



same form, the only difference being that A in the first 
equation is replaced by A in the second, and by A in the 
third one. In what follows we are going to use the symbol 
A, understanding that 

i) for 7-unstable nuclei it is given by A = r(r + 3), 

ii) for axially symmetric prolate d efor med nuclei it 
should be replaced by A, given in Eq. |43| , and 

iii) for triaxial nuclei it should be replaced by A, given 
in Eq. 



Eq. 

tions 



52) 



32 ) can be simplified by performing the deriva- 



2£ 





U'Y 

8 



4 



f£ 

/3 



with 



■^pA£ + e£ - v e ff£ = 0, 



(53) 



1 (I 



2 V 2 



-6 



1 



- A (/ ) 



(54) 



The difference in the numerical coefficient of /' observed 
in comparison to Eq. (2.27) of Ref. [9] is due to the 
different dimensionality of the space used in each case. 
Setting 



R(P) 



Eq. ( 53 ) is put into the form 
d 



HR 

where 



v/7 



l eff 



dp 



,J eff 



2 

/I R 



2u eff R = 2eR, 



P + Pff , P 



P* 



2p 2 



A. 



(55) 



(56) 



(57) 
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V. THE DAVIDSON POTENTIAL 

Up to now no assumption about the specific form of 
the potential u(f3) and the deformation function f(/3) has 
been made. We are now going to consider the special case 
of the Davidson potential [13] 



u{p) = p 2 



Pt 



(58) 



where the parameter /3q indicates the position of the min- 
imum of the potential. The special case of Pq = corre- 
sponds to the simple harmonic oscillator. 

Based on the results for the 3-dimensional harmonic 
oscillator reported in Ref. [T^], we are also going to con- 
sider for the deformation function the special form 



/(/?) = 1 + ap 2 



a< 1. 



(59) 



This choice is made in order to lead to an exact solution. 
Its physical implications will be discussed in Section 11. 
Using these forms for the potential and the deforma- 



tion function in Eq. (57) one obtains 



2 . , k -i 



2u eS = ki0 z + fc + 



/3 2 ' 



(60) 



where 

ki = 2 + a 2 [5(l - S - A) + (1 - 2<5)(1 - 2A) + 6 + A], 
fc = a[5(l-tf-A) + 8 + 2A], 
jfe_i = 2 + A + 2/3^. 

(61) 



VI. DEFORMED SHAPE INVARIANCE 

Our task now is to find the eigenvalues and eigen- 
functions of the Hamiltonian of Eq. ( [56] ). This can be 
achieved by imposing shape invariance j!2j. which is an 
integrability condition guaranteeing that exact solutions 
of the Hamiltonian of Eq. (56 1 can be found. The use of 



shape invariance in the framework of SUSYQM is equiv- 
alent pj] to the well known factorization method of the 
Schrodinger equation, introduced 60 years ago by Infeld 
and Hull |25j . In other words, we are now going to use a 
mathematical technique allowing us to find the solutions 
ofEq. (1561). 



In its simplest form in an one-dimensional space, shape 
invariance can be described as follows [33] ■ Two poten- 
tials Vi and V 2l which are supersymmetric partners, are 
in general different functions of x. They are called shape 
invariant if they satisfy the condition 



V 2 (x;ai) = Vi(x;a 2 ) + R(ai), 



(62) 



where <zi, a 2 are sets of parameters independent of x, 
with a 2 being a function of a±, and the remainder R(ai) 
is also independent of x. In other words, the two po- 
tentials have the same functional dependence on x, the 



difference being in the values of the parameters appear- 
ing in each of them, and in their relative displacement by 
the remainder i?(ai). Furthermore, it is known that the 
shape invariance condition of Eq. ( 62 ) can be written in 
the operator form 



A( ai )A^( ai ) = A\a 2 )A(a 2 ) + R( ai ), 



(63) 



where A and A^ are the operators corresponding to the 
supersymmetric partners Hi = A^A and H 2 = AA\ 
Solving the Schrodinger equation for H\ by this method, 
one obtains as a "bonus" the solution of H 2 as well. 

In the present case, the concept of shape invariance 
has to be generalized, as described in detail in Ref. |12) . 
since the mass depends on the deformation, resulting in 
a deformed shape invariance condition. Instead of two 
Hamiltonians, one has a series of many Hamiltonians. 
We are interested in solving the Schrodinger equation for 
the first of them, which will be Eq. ( 56 1 . 



H in Eq. ( 56 ) may be considered as the first member 
Hq = H of a hierarchy of Hamiltonians 



Hi = A+Ar + J2^, * = 0,1,2... 

3=0 



where the first-order operators [12] 



(64) 



Af =A ± { l ii,Vi) = T^]j- p ^] + W{^, Vt] P) (65) 



satisfy a deformed shape invariance condition 



A7A1 



A i+1 A i+1 + Ei+i, 



* = 0,1,2, 



(66) 



with Ei, i = 0, 1, 2, . . . , denoting some constants. (Note 
that the parameters A and 11 of [12] have been changed 
into \x and v, respectively.) 

In other words, the superpotential W(fi, v\ (3) fulfils the 
two conditions 

W 2 (p, v\ P) - f(fi)W'Qi, v- /3) + £ = 2 Ueff (/3) (67) 

and 

IU 2 ( Ml ,^ ;/?) + / (,3) IU'( Ml ,^;/3) 
= W 2 (jm +1 , u i+1 ; p) - f{fi) W'(f J L i+1 , v l+1 ■ p) + E i+1 , 

i = 0,1,2,..., (68) 

where /io = fx, vo = v, and a prime denotes derivative 
with respect to p. 



In the case of the effective potential given in Eq. ( 60 1 , 
W(}i, v; P) is a class 2 superpotential 



W(n,v t p) = tuf>(P) + 



1 



<kpv 



MP) = a. 



(69) 
(70) 



which means that Eqs. (3.9) and (3.10) of [H] read 

1 A 
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{A' /p 2 ) + B< 



(71) 



with A = -1, B = 0, A' = 0, and B' = -a. 



Inserting Eqs. (69) and (70) in (67), we obtain 



h(3 2 + k Q + 

P 



2 



(72) 



which is equivalent to the three equations 

/i(/i+ 1) = fc-i, v{y — a) = k\, 2\iv + \ia — v + Eq = fco- 

(73) 

Their solutions read 

/*=^(-l±Ai), i/ = |(l±A 2 ), 



Ai = v/l + 4fc-i, 



A 2 = Wl + 4-i, 



(74) 



provided 1 + 4fci/a 2 > (note that 1 + 4fc_i is always 
positive). As we shall show in Sec. VIII A the conditions 
ensuring that the ground-state wavefunction is physically 
acceptable select the lower sign for /z and the upper one 
for v: 



1 



1 + AO 



v=~(l + A 2 ). (75) 



Inserting next Eqs. ( 69 ) and ( 70 1 in Eq. ( 68 ) , we get 

2 



M*+l i a\ n , q2\ ( M«+l 
+ u i+ i/3 - (1 + a/3 ) 



leading to the three conditions 



P 2 



(76) 



Mi (Mi ~ 1) = Mi+i(Mi+i + 1)) 
i/,-(t/j + a) = - a), 

2jLt»i/i - fiia + V{ = 2(i i+1 L> i+1 + (ii+ia - v i+x + e l+ i. 

(77) 



Their solutions are 

Mi+i = Mi - 1) = "i + fl i 

and 



(78) 



Note that there are other solutions for Hi+i and 
namely = — /x, and t^ + i = — z/£, but the alternating 
signs would not be compatible with physically acceptable 
excited-state wavefunctions. Finally, the iteration of ( 78 ) 
leads to 



Vi = V + Ml. 



(80) 



VII. ENERGY SPECTRUM 



The energy spectrum of Eq. ( 56 1 is therefore given by 

n 

£ " = 2 £i 



i=0 



fco — 2fl n l/ n 



2 ^ (J,i + lln + 2 ^ Vj + l/„ 



i=Q 



i=0 



= I [fco ^ 2ii^ — azx + v — 4(a/x — v)n + Aan 2 ] . 



(81) 



On taking (75) into account, this can be rewritten as 

(82) 



e„ = \ [fco + §a(3 + 2Ai + 2A 2 + A 1 A 2 



2a(2 + Ai + A 2 )rx + 4arT 



n = 0,1,2,. 



Equation ( 82 ) only provides a formal solution to the 



bound-state energy spectrum. The range of n values is 
actually determined by the existence of corresponding 
physically acceptable wavefunctions. The relevant con- 
ditions will be considered in the next section. 

We quote here the final results for the spectra, which 
will be used for comparison to experiment. One has 



ro 
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a+-(l-S- X)a+ - yja 2 + 4fci 



+ 4fc_i + - + 4fci)(l + 4A_i) + aA, 



(83) 



ei = (a 



+ 4a+ \/a 2 +4fci + ay/l + 4fc-i, (84) 



e 2 = e + 12a + 2^/a 2 + 4fci + 2a^l + 4fc_ l5 (85) 



where fci, fc_i are given by Eq. (61 1, in which A has the 
form explained in subsec. |IVD| 



The ground state band is obtained from Eq. ( 83 1 , while 



the quasi- ft band is obtained from Eq. (84), and the 
quasi-/3 2 band is obtained from Eq. (85). 



In the special case of a = (no dependence of the mass 
on the deformation) one easily obtains 



ei = eo 



2^2, 



£2 = eo 



4^2, 



(86) 



= 2(//ji/j — /Xj_|_ii/j_|_i) — (/i,-+/i,_|_i)a+i'j+J/j-|-i. (79) i.e. the /3-bandheads become equidistant. 
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VIII. WAVE FUNCTIONS 

To be physically acceptable, the bound-state wavefunc- 
tions should satisfy two conditions [12]: 
(i) As in conventional (constant-mass) quantum mechan- 
ics, they should be square integrable on the interval of 
definition of u e s, i.e., 



(61 1 will be positive and A2 in (74 1 will be greater than 
1. For sufficiently large values of a, however, both sign 
choices might be acceptable. Since among two acceptable 
wavefunctions, it is customary in quantum mechanics to 
choose the most regular one (see, e.g., [33] and references 
quoted therein), we assume the upper sign for v, thus 
getting Eq. (75]). 



d(3\R n {t3)\ 2 <oo. 



(87) 



(ii) Furthermore, they should ensure the Hermiticity of 
H. For such a purpose, it is enough to impose that the 
operator \ff {d/ dj3)\ff be Hermitian, which amounts to 
the restriction 



\R n ((3)\ 2 f(f3) -> 
or, equivalently, 



for fi — > and (3 — > 00, 



\Rn(/3)\ 2 -> for/3^0 
\Rn{P)\ 2 P 2 ^0 for/3 



and 
> 00. 



(89) 



As condition ( 89 1 is more stringent than condition ( 87 1 , 



we should only be concerned with the former. 



A. Ground-state wavefunction 

The ground-state wavefunction, which is annihilated 
by A~, is given by Eq. (2.29) of [12 as 



Nn 



: exp 



" W ^di\, (90) 



where N is some normalization coefficient. Here 



P l + a/3 2 J 



1 



= /iln/3+ — (^-^a)ln(l + a^ 2 ). (91) 
2a 



Hence 



(92) 



For p — > 0, the function \R (P)\ 2 behaves as (3~ 2 ^. 
Condition ( [89] ) imposes that —2/j, > or fi < 0. Since 
defined in Eq. (61 1, is greater than 2, it follows that 



Ai, defined in (74), is greater than 3, so that the upper 



sign choice for \x in (74) would lead to fx > 1. As this is 
not acceptable, we have to take the lower sign for which 
H < -2. 

Jbr (3 -> 00, \R (P)\ 2 P 2 behaves as P~ 2v / a . Condition 



(89) therefore imposes that v > 0. This restriction is 



surely satisfied by the upper sign choice for v in ( 74 1 . 
For the lower one, it is not fulfilled if we restrict our- 
selves to small enough values of a because then k\ in 



B. Excited-state wavefunctions 

According to Eqs. (2.30), (3.20) and (3.21) of 12J, the 
excited-state wavefunctions are given by 

R„{/3) = R n (n, v\ 0) (x f3~ n R (p n , v n ;P)P n (n, v; y), 

(93) 

where P n (fj,,v;y) is an nth-degree polynomial in y, sat- 
isfying the equation 



P n+1 (fi,v;y) = -2y(l + ay)—P n ((i 1 ,v 1 ;y) 
+ \p n +i + " + n + (u n+ i + v + na)y]P n (fx-L,i>i;y), 



(94) 



with the starting value Pq((J., v\ y) = 1. 
From Eqs. (flot) and pfb, it follows that 



KR (^v;P)P n r n , 



so that Eq. ( 93 ) becomes 



R n (/3) cx R (/3)f- n P n (fi,v;y). 



(95) 



(96) 



It is then clear that R n {P) satisfies condition |89[) for any 
n = 1, 2, . . . , since Ro(ft) does. 



It now remains to solve Eq. (94). For such a purpose, 



let us make the changes of variable and of function 



t = 1 



-1 + ay 



f 1 + ay ' 
Pni^vw) =C n f n Q n (pL,v;t), 



(97) 



where C n is some constant. From definition (97), it fol- 
lows that Q n ([i,i>;t) an nth-degree polynomial in t. We 
successively get 



1 + t 1 2 
a(l — t) 1 — t 



d a 2 d 

dy ~ 2^ ~ 1 ' dt' 



(98) 

It is then straightforward to show that Eq. ( 94 ) becomes 

dt (99) 





V 




-.).] 




h 






a 
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On taking into account that the Jacobi polynomi- 
als satisfy the backward shift operator relation (see 
Eq. (1.8.7) of [35]) 



2(n+l)P&£> (*) = j-(l-z 2 ) ( 



d 

' dx (100) 

+ [a -13 + {a + /3 + 2)x}} P 7 [ a+hl3+1 \x), 

we see that Q n (/i, v; t) is actually some Jacobi polynomial 



1 ( At)I /;i)=pi" *' " h \t) = pi %2 - %1 \t), (101) 



provided we choose 



%±i = 2(n + l), G = l, 



(102) 



or, in other words, C n — 2 n n\. 

We therefore conclude that the wavefunctions are given 
by 



iVo 



j—{v—ixa+a)/(2a)p\ a i< V \ ) ^\ 



(103) 



i?„(/3) = iVnjS (l+A 1 )/2 / -l-(A 1 +A 2 )/4p j (A 2 /2,A 1 /2) (f)) 

-l + a/3 2 



t = 



1 + a^ 2 ' 



(104) 



where N n is some normalization coefficient. 

The Jacobi polynomials appearing in the wave func- 
tions of the ground state band (n = 0), the quasi- f3\ 
band (n = 1), and the quasi-/?2 band (n = 2), needed for 
the calculation of the relevant B(E2) transitions, read 



Pf> f> \x) = l, 



(105) 



P^\x) = |[2(a+l) + (a + ^ + 2)(x-l)], (106) 



(107) 



P^\x) = ~[4(a + l)(a + 2) 

+ A{a + (3 + + 2){x - 1) 
+ (« + /? + 3) (a + 13 + 4)(x - l) 2 ]. 

IX. NORMALIZATION COEFFICIENT 



To calculate iV„, let us first express the whole wave- 
function R n in terms of t: 

Rn = iV„Z/ (1+Al)/4 (l + ay)- 1 -^+ A ^/ 4 P 7 [ A ^ 2 - A ^ 2 \t) 
= iV n 2 _1_(Al+A2)/4 a _(1+Al)/4 (l + t) {1+Al),A 
(1 - t)( 3+A2)/4 P {A2/2 ' Al/2) (t). 

(108) 



Now, on taking into account that 
dy dt 



df3 



we obtain 



2^y v^(l-t)3/3(l + t)i/2 



(109) 



|i?„| 2 d/3 = \N n \ 2 2- 2 -^ +A ^/ 2 a- 1 - A ^ 2 



+ i 



(l-t) A2/2 (l + t) 



Ax/2 



P, 



(A 2 /2,A 1 /2) 



(*) 



(110) 



in terms of the normalization integral of Jacobi polyno- 
mials [56] , 

Hence the normalization condition reads 
|iV n | 2 2- 2 -( Al+A2 ^ 2 a - 1 - Al / 2 

2 (A 1+ A 2 )/2+l r ( n 

(2n + + 1) nir ( 

and leads to 



+ ^ + l)T(n+ ^ + 1) _ (in) 
l)n!rfn + ^±Aa + l) ' 



A,, = (^2a Al ^ 2+1 (2n+^±^ + l 



1/2 



r (r 



A t + A 2 



1 



1/2 



(112) 



r(n + ^- + i)r(- 1 A 



A way of avoiding numerical problems when having to 
handle T(x) functions with large x is given in Appendix 
1. 



X. B(E2) TRANSITION RATES 

B(E2) transition rates 

5 \(g'L'\\T( E ^\\gL)\ 2 



B{E2\qL -> g'L 



16tt 2L + 1 



(113) 



where g stands for quantum numbers other than the 
angular momentum L, can be calculated using the 
quadrupole operator T^ 2 ) and the Wigner-Eckart theo- 
rem in the form 



(g'L'M'\T^\gLM) 

A_ (L2L'\MnM') (q'L' \\T^ \\qL). 



(114) 
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A. B(E2)s for 7-unstable nuclei 

The calculation is carried out exactly as in Ref. |37) . 
using the quadrupole operator [31] 



T (E2) = A/3 



P^(^)co S7 



+ ^(^ 2 2 (^)+^ 2 (^))sin 7 



(115) 
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where A is a scale factor. 

The results of Ref. [37] need not be repeated here. The 
only difference is that in the radial integral (see Eq. (21) 
of Ref. [57]) the wave functions i?„ jT (/3) appear 



n' ,r+l;n,r 



(116) 

v {l3)R n , T {P)df3. 



The r dependence of the wave functions R n (/3) of Eq. 
(104 1 is contained in Ai, A2, known from Eq. (74) to 



contain k\, k—\, which in turn are known from Eq. (|61|) 
to contain A = r(r + 3). 



B(E2)s for axially symmetric prolate deformed 
nuclei 



The quadrupole operator is again given by Eq. (1151. 
The calculation is carried out exactly as in Ref. |18|. the 
results of which need not be repeated here. The only 
difference is that in the radial integral (see Eq. (B5) of 
Ref. [18]) the wave functions R n ,L,{P) appear 



(117) 



B n ,L,n',L' 

(3R n , L (/3)Rn>,L>(P)df3 



The L dependence of the wave functions R n {@) of Eq. 
(104 1 is contained in Ai, A2, known from Eq. (74) to 



contain ki, k—\, whi ch i n turn are known from Eq. (61) 
to contain A of Eq. ( 43 ) . 



C. B(E2)s for triaxial nuclei with 7 = n/6 

The calculation is carried out exactly as in Ref. 
using the quadrupole operator 



Tl E 2 ) = A[3 



V 



(2) in 



) cos 7 — 



2tt 

T 



(118) 



where A is a scale factor, while the quantity 7 — 27r/3 in 
the trigonometric functions is obtained from 7 — 2irk/3 
for k = 1, since in the present case the projection a along 
the body-fixed i'-axis is used. 

The results of Ref. [35] need not be repeated here. The 
only difference is that in the radial integral (see Eq. (14) 
of Ref. [38]) the wave functions R n , a ,L{fl) appear 



Ip(n, L, a; n , L\ a') 



Kn, a .Lmn>,L>, a -{P)P A dP 



PRn, a ,L{P)Rn',a>MPW- 



(119) 



The a, L dependence of the wave functions R n (f3) of Eq. 
(104) is contained in Ai, A2, known from Eq. (74) to 
contain k±, whi ch i n turn are known from Eq 
to contain A of Eq. (51 ). 



(61) 
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FIG. 1: The function /3 2 /f 2 (/3) = /3 2 /(l + afi 2 ) 2 , to which mo- 
ments of inertia are proportional as seen from Eq. |9|, plotted 
as a function of the nuclear deformation j3 for different values 
of the parameter a. See Section XI for further discussion. 
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FIG. 2: Moments of inertia 0(1/) for the ground state band, 
given by Eq. ( |120[ ) and normalized to 0(2), are shown for the 
specific values of /3o = 2 and c = 5, and varying parameter a. 
See Section XI for further discussion. 



XI. NUMERICAL RESULTS 

From Eq. ^ it is clear that in the present 
case the moments of inertia are not proportional to 
/3 2 sin 2 ( 7 -27r/c/3) but to (/3 2 // 2 (/3)) sin 2 (7 - 2tt/c/3). 
The function /3 2 // 2 (/3) is shown in Fig. 1 for different 



12 



values of the parameter a. It is clear that the increase 
of the moment of inertia is slowed down by the function 
/(/?), as it is expected as nuclear deformation sets in [3]. 

The effect of the deformation-dependent mass on the 
moments of inertia can be seen in Fig. 2, where the mo- 
ments of inertia [3] for the ground state band 



G(L) 



2L- 1 



E{L) -E{L-2) 



(120) 



normalized to 6(2), are shown in the case of axially sym- 
metric prolate deformed nuclei, for the specific values of 
/3 = 2 and c = 5, and varying parameter a. It is clear 
that the rapid increase of the moments of inertia with L, 
seen for a = 0, is gradually moderated by increasing a. 



A. Spectra of 7-unstable nuclei 

Rms fits of spectra have been performed, using the 
quality measure 



Tr=i(^(e^)-^(^)) 2 



(n-l)E(2+ 



(121) 



The theoretical predictions for the levels of the ground 
state band are obtained from Eq. ( 83 ) , while the levels of 



the quasi- /?i band are obtained from Eq. (84 1. The levels 



of the quasi-71 band are obtained through their degen- 
eracies to members of the ground state band, mentioned 
below Eq. 

The results shown in Table 1 have been obtained for 
S = A = 0. (The Xe and Ba isotopes have already been 
considered in Ref . [23] ■ ) One can easily verify that differ- 
ent choices for 5 and A lead to a renormalization of the 
parameter values a and /3o, the predicted energy levels 
remaining exactly the same. 

Concerning the physical content of the parameter a, 
it is instructive to consider in detail in Table 1 the 
Xe isotopes (already discussed in Ref. [23] ), known 
|39j to lie in a 7-unstable region. They extend from 
the borders of the neutron shell ( 134 Xe§o is just below 
the N=82 shell closure) to the midshell ( 120 Xe 66 ) and 
even beyond, exhibiting increasing collectivity (increas- 
ing Ri/2 — E(4i)/E(2f) ratios) from the border to the 
mishell. Moving from the border of the neutron shell to 
the midshell, the following remarks apply 

i) 134 Xe and 132 Xe are almost pure vibrators. There- 
fore no need for deformation dependence of the mass ex- 
ists, the least square fitting leading to a = 0. Further- 
more, no /3q term is needed in the potential, the fitting 
therefore leading to fig = 0, i.e., to pure harmonic be- 
haviour. 

ii) In the next two isotopes ( 130 Xe and 128 Xe) the 
need to depart from the pure harmonic oscillator becomes 
clear, the fitting leading therefore to nonzero fto values. 
However, there is still no need of dependence of the mass 
on the deformation, the fitting still leading to a = 0. 



iii) Beyond 126 Xe both the /?o term in the potential and 
the deformation dependence of the mass become neces- 
sary, leading to nonzero values of both /3 and a. 

Other chains of isotopes also show similar behavior. 



B. Spectra of axially symmetric deformed nuclei 

Fits of spectra of deformed rare earth and actinide nu- 
clei are shown in Table 2. The energy levels of the ground 

and 



state band and the (3\ band (both having n 7 
K = 0) are obtained from Eqs. ([83]) and (84| respec- 



tively, while the levels of the 71 band are obtained from 



Eq. (83) with n 7 = 1 and K = 2. Again, the choice 



5 = A = has been made, and it is seen that different 
choices for S and A lead to a renormalization of the pa- 
rameter values a, /3o, and c, the predicted energy levels 
remaining exactly the same. 

The quality of the fits obtained can also be seen in 
Table 3, where the calculated energy levels of 162 Dy and 
238 U are compared to experiment. 

The following remarks apply. 

1) Both the bandheads and the spacings within bands 
are in general well reproduced. This is particularly true 
for the ground state and the 71 bands. The deviation 
in the gsb of 162 Dy reaches 0.6% at L = 18, while in 
the gsb of 238 U it reaches 1.7% at L = 30. The exper- 
imental levels of the 71 band of 162 Dy (up to L = 14) 
extend over 28.4 energy units, while the corresponding 
theoretical predictions spread over 28.7 units, the dif- 
ference being of the order of 1%. Similarly in 238 U the 
experimental spread of the 71 band (up to L = 27) is 
89.1 energy units, while the theoretical one is 87.3 units, 
the difference being of the order of 2%. 

2) However we remark that the theoretical level spac- 
ings within the f3\ bands are larger than the experimen- 
tal ones. This should be attributed to the shape of the 
Davidson potential, which raises to infinity at large /?, 
pushing f3 bands higher and increasing their interlevel 
spacing. It is known that this problem can be avoided 
by using a potential going to some finite value at large ft 
[31], like the Morse potential [32] • 



C. B(E2)s of 7-unstable nuclei 

B(E2)s within the ground state band, as well as in- 
terband B(E2)s for which experimental data exist for 
several nuclei, have been calculated using the procedure 
described in subsec. X.A . The results are shown in Table 
4, the overall agreement being good. 



D. B(E2)s of axially symmetric deformed nuclei 

B(E2)s within the ground state band, as well as inter- 
band B(E2)s for which experimental data exist for sev- 
eral nuclei, have been calculated using the procedure de- 
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scribed in subsec. X.B . The results are shown in Table 5. 
The overall agreement is good for transitions within the 
ground state band (gsb), as well as for transitions con- 
necting the 71 band to the gsb, while transitions from 
the 0\ band to the gsb tend to be overpredicted. One 
should remember at this point that the 0\ band was the 
one giving poor results also in the case of the spectra, in 
subsec. XI.B . 



XII. CONNECTION TO EARLIER WORK 

It is instructive to examine the relation between the 
present approach and earlier numerical work. 

1) The formalism of subsection III.B clarifies the re- 
lation between the present approach and the numerical 
solution of Kumar and Baranger |26) , who used a matrix 
of the form (14) with gij, i,j = 1,2,3 the same as in Eq. 
( 15 1, but with 



#44 — Bqq, g55 = B 2 '2', 945 — .954 — #02' i (122) 

where Boo, B^i', B02', as well as the moments of inertia 
Ji (i = 1, 2, 3) and the potential V have been treated as 
seven arbitrary functions of the variables 0q — 0cosj 
and 02' — sin 7 [while in the Bohr formulation [1 
ao = /3cos7 and 02 = 0siwy/y/2 are used]. On one 
hand, the present solution is a special case of Ref. , 
since it contains no non-diagonal terms c/45 = (754. On the 
other hand, in the present approach the above mentioned 
quantities are interrelated by the overall symmetry in a 
specific way, greatly reducing the number of free parame- 
ters (down to two or three in total). It should be pointed 
out that the functional dependence of the mass on the de- 
formation for the potential used is dictated by SUSYQM. 
Therefore, the successful prediction of the behavior of the 
moments of inertia, for example, provides credit for the 
present approach. What we see, independently of the 
parameter values, is that the increase of the moments of 
inertia as a function of deformation is moderated by the 
/ 2 factor, which can be seen as a result of the dependence 
of the mass on the deformation, or, alternatively, as seen 
in subsection III.B, as a result of using a curved space. 

2) It should be pointed out that in Ref. [9] the equiva- 
lence between the position dependent mass case and the 
curved space approach has been established in the special 
case of k — 2 and S = A = (see Eq. ^ for the mean- 
ing of the symbols), which represents the BenDaniel and 
Duke Hamiltonian 143 



XIII. CONCLUSION 

In the present work analytical solutions are obtained 
for a Bohr Hamiltonian in which the mass has been al- 
lowed to depend on the deformation. 

From the mathematical point of view, this is achieved 
through the use of techniques of supersymmetric quan- 
tum mechanics [101 111] , involving a deformed shape in- 
variance condition [12 . Analytical expressions for the 
spectra and wave functions have been obtained. 

From the physics point of view, spectra and B(E2) 
transition rates have been calculated for 7-unstable, axi- 
ally symmetric prolate deformed, and triaxial nuclei, im- 
plementing the usual approximations in each case, and 
compared to experimental data for the first two cases. 
The main new result is that the dependence of the mass 
on the deformation moderates the increase of the moment 
of inertia with the deformation, removing an important 
drawback j3] of the model. It should be emphasized that 
the functional dependence of the mass on the deforma- 
tion for the potential used is dictated by SUSYQM, thus 
the correction in the behavior of the moments of inertia 
is a general effect, independent of any specific parameter 
value combinations. 

However, certain discrepancies with experimental data 
remain, especially related to the /?i-band and its inter- 
band transitions. It should be remembered at this point 
that in the present study separation of variables has been 
achieved by assuming that the potential either is indepen- 
dent of the 7- variable, or it has the exactly separable form 
of Eq. (34 1. Furthermore, the approximations related to 
Eqs. (35) and (44) have been implemented. Recently, the 



numerical solution of the Bohr Hamiltonian for any value 
of and 7, avoiding all these approximations, has been 
achieved in the framework of the powerful algebraic col- 
lective model [4"4T(4l)] . The detailed study of discrepancies 
from experimental data both in the SUSYQM framework 
and in the context of the algebraic model, especially for 
multi-phonon excitations [57], could shed light on the 
origins of these discrepancies. 

As it has already been mentioned, the form of the 
dependence of the mass on the deformation is dictated 
by SUSYQM for the potential used in the degree 
of freedom. In the present work, the Davidson poten- 
tial has been used, called the Deformation Dependent 
Mass (DDM) Davidson model. The application of the 
SUSYQM approach to the Bohr Hamiltonian with the 
Kratzer potential [48, 49 is receiving attention. 



H 



BD 



-V/ 2 V + U SD . 



(123) 



This resembles the collective Hamiltonian 

H ^ = -Y^-i qi [M(q)i ^ 1 ^ +v(q) (124) 

used by Libert et al. [27] in mean field calculations, in 
which a tensor mass appears. 
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Appendix 1 



r n 



A, 



r(J + (n + r + l)) 



When using Eq. (112) in numerical calculations, prob- 



(130) 



lems can appear because of T(x) functions with large x. 
These problems can be avoided by using Eq. 6.1.16 of 
Ref. [II] 

T{n + z) = (n-l+z)(n-2 + z) . . . (l + z)T(l + z). (125) 
In the normalization factors we need the ratio of 



over 



Ai + A 2 

r n H h 1 



A2 

2 



(126) 



(127) 



= (I-l + n + r + I) (I -2 + n + r + l). 
(1 + n + r + l)r(l + n + r + 1). 

Their ratio becomes 



r(n+^|^ + l) = {I- l + n + Aj/2 + r + l) 

r(n+^ + l) ~ (J-l + n + r+1) 
(/-2 + n + Ai/2 + r + l) (1 + n + Ai/2 + r + 1) 

(7 -2 + n + r + l) "' (1 + n + r + l) 
r(l + n + Ai/2 + r + l) 



r(l + n + r + l) 
Ar/2 



Let us call I the integer part of A 2 /2 and r the rest of 
it, i.e., 



(I -1 + n + r + l) 
Ai/2 



Ai/2 



= I + r. 



(128) 



Then we have 
r n H 



= r(J + (n + Ai/2 + r + l)) 
= (J-l + n + Ai/2 + r + l) 

(I - 2 + n + Ai/2 + r + 1) . . . (1 + n + Ai/2 + r + 1) 
r(l + n + Ai/2 + r + l), 



(7 -2 + n + r + l)/ \ (1 + n + r+l), 
r(l + n + Aj/2 + r + l) 
T(l + n + r + l) 

(131) 



in which one does not have to calculate T(x) functions 
with large x. The only large numbers appear in denom- 
inators of fractions accompanying 1, which do not pose 
any problem. 
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TABLE I: Comparison of theoretical predictions of the 7-unstable Bohr Hamiltonian with /3-dependent mass (with S — A = 0) 
to experimental data [40] of rare earth and actinide nuclei with R4/2 < 2.6 and known 0^" and 2+ states. The R. 



M/2 = 



E(4f)/E(2f) ratios, as well as the P and 7 bandheads, normalized to the 2f state and labelled by R0/2 = E(0p)/ E(2f ) and 
R2/2 = I E{2^ ) respectively, are shown. Po and a are free parameters, related to the Davidson potential [Eq. (58)] and 

to the dependence of the mass on the deformation [Eq. (59 1]. The angular momenta of the highest levels of the ground state, 



P and 7 bands included in the rms fit are labelled by L g 
involved in the fit and a is the quality measure of Eq. 
mentioned below Eq. dl21 
discussion. 



Lp, a nd L 7 respectively, while n indicates the total number of levels 
(121 1. The theoretical predictions are obtained from the formulae 
The Xe and Ba isotopes have already been considered in Ref. [53]. See subsec. |XI A| for further 



nucleus R^n 
exp 


R 


1/2 
th 


Ro/2 

exp 


Rq/2 

th 


R2/2 
exp 


R2/2 
th 


Po 


a 


L 9 


Lp 




n 


a 


98 Ru 


2 


14 


2 


14 


2.0 


2.4 


2.2 


2 


1 


0.99 


0.020 


24 





4 


15 


0.277 


100 Ru 


2 


27 


2 


24 


2.1 


2.7 


2.5 


2 


2 


1.19 


0.048 


28 





4 


17 


0.315 


102 Ru 


2 


33 


2 


20 


2.0 


2.4 


2.3 


2 


2 


1.05 


0.059 


16 





5 


12 


0.364 


104 Ru 


2 


48 


2 


34 


2.8 


3.0 


2.5 


2 


3 


1.40 


0.083 


8 


2 


8 


12 


0.429 


102p d 


2 


29 


2 


24 


2.9 


2.3 


2.8 


2 


2 


1.08 


0.081 


26 


4 


4 


18 


0.326 


104p d 


2 


38 


2 


21 


2.4 


2.6 


2.4 


2 


2 


1.15 


0.034 


18 


2 


4 


13 


0.397 


106p d 


2 


40 


2 


16 


2.2 


2.2 


2.2 


2 


2 


0.91 


0.062 


16 


4 


5 


14 


0.409 


108p d 


2 


42 


2 


26 


2.4 


2.3 


2.1 


2 


3 


1.09 


0.103 


14 


4 


4 


12 


0.318 


uo pd 


2 


46 


2 


31 


2.5 


2.0 


2.2 


2 


3 


0.99 


0.195 


12 


10 


4 


14 


0.354 


112p d 


2 


53 


2 


29 


2.6 


2.5 


2.1 


2 


3 


1.21 


0.086 


6 





3 


5 


0.485 


114p d 


2 


56 


2 


31 


2.6 


2.8 


2.1 


2 


3 


1.30 


0.076 


16 





11 


18 


0.722 


116p d 


2 


58 


2 


36 


3.3 


3.4 


2.2 


2 


4 


1.52 


0.062 


16 





9 


16 


0.609 


106 Cd 


2 


36 


2 


25 


2.8 


2.9 


2.7 


2 


3 


1.28 


0.028 


12 





2 


7 


0.268 


108 Cd 


2 


38 


2 


14 


2.7 


2.2 


2.5 


2 


1 


0.91 


0.041 


24 





5 


16 


0.528 


110 Cd 


2 


35 


2 


08 


2.2 


1.9 


2.2 


2 


1 


0.00 


0.061 


16 


6 


5 


15 


0.415 


112 Cd 


2 


29 


2 


05 


2.0 


1.9 


2.1 


2 





0.00 


0.033 


12 


8 


11 


20 


0.523 


114 Cd 


2 


30 


2 


06 


2.0 


1.9 


2.2 


2 


1 


0.00 


0.041 


14 


4 


3 


11 


0.418 


116 Cd 


2 


38 


2 


16 


2.5 


2.7 


2.4 


2 


2 


1.14 


0.000 


14 


2 


3 


10 


0.387 


118 Cd 


2 


39 


2 


19 


2.6 


2.9 


2.6 


2 


2 


1.21 


0.002 


14 





3 


9 


0.429 


120 Cd 


2 


38 


2 


20 


2.7 


2.9 


2.6 


2 


2 


1.22 


0.006 


16 





2 


9 


0.412 


118 Xe 


2 


40 


2 


32 


2.5 


2.6 


2.8 


2 


3 


1.27 


0.103 


16 


4 


10 


19 


0.319 


120 Xe 


2 


47 


2 


36 


2.8 


3.4 


2.7 


2 


4 


1.51 


0.063 


26 


4 


9 


23 


0.524 


122 Xe 


2 


50 


2 


40 


3.5 


3.3 


2.5 


2 


4 


1.57 


0.096 


16 





9 


16 


0.638 


124 Xe 


2 


48 


2 


36 


3.6 


3.5 


2.4 


2 


4 


1.55 


0.051 


20 


2 


11 


21 


0.554 


126 Xe 


2 


42 


2 


33 


3.4 


3.1 


2.3 


2 


3 


1.42 


0.064 


12 


4 


9 


16 


0.584 


128 Xe 


2 


33 


2 


27 


3.6 


3.5 


2.2 


2 


3 


1.42 


0.000 


10 


2 


7 


12 


0.431 


130 Xe 


2 


25 


2 


21 


3.3 


3.1 


2.1 


2 


2 


1.27 


0.000 


14 





5 


11 


0.347 


132 Xe 


2 


16 


2 


00 


2.8 


2.0 


1.9 


2 





0.00 


0.000 


6 





5 


7 


0.467 


3 Xe 


2 


04 


2 


00 


1.9 


2.0 


1.9 


2 





0.00 


0.000 


6 





5 


7 


0.685 


130 Ba 


2 


52 


2 


42 


3.3 


3.2 


2.5 


2 


4 


1.60 


0.118 


12 





6 


11 


0.352 


132 Ba 


2 


43 


2 


29 


3.2 


2.8 


2.2 


2 


3 


1.29 


0.059 


14 





8 


14 


0.619 


134 Ba 


2 


32 


2 


16 


2.9 


2.7 


1.9 


2 


2 


1.12 


0.000 


8 





4 


7 


0.332 


136 Ba 


2 


28 


2 


00 


1.9 


2.0 


1.9 


2 





0.00 


0.000 


6 





2 


4 


0.250 


142 Ba 


2 


32 


2 


38 


4.3 


4.3 


1.0 


2 


4 


1.72 


0.028 


14 





2 


8 


0.609 


134 Ce 


2 


56 


2 


34 


3.7 


3.9 


2.4 


2 


3 


1.59 


0.019 


34 


2 


8 


25 


0.527 


136 Ce 


2 


38 


2 


11 


1.9 


2.1 


2.0 


2 


1 


0.82 


0.034 


16 





3 


10 


0.457 


138 Ce 


2 


32 


2 


00 


1.9 


2.0 


1.9 


2 





0.00 


0.000 


14 





2 


8 


0.314 


140 Nd 


2 


33 


2 


05 


1.8 


1.9 


1.9 


2 


1 


0.00 


0.037 


6 





2 


4 


0.192 


148 Nd 


2 


49 


2 


36 


3.0 


2.8 


4.1 


2 


4 


1.38 


0.110 


12 


8 


4 


13 


0.764 


140 Sm 


2 


35 


2 


29 


1.9 


1.9 


2.7 


2 


3 


0.92 


0.196 


8 





2 


5 


0.207 


142 Sm 


2 


33 


2 


06 


1.9 


1.9 


2.2 


2 


1 


0.33 


0.044 


8 





2 


5 


0.147 


142 Gd 


2 


35 


2 


21 


2.7 


2.8 


1.9 


2 


2 


1.20 


0.020 


16 





2 


9 


0.231 


144 Gd 


2 


35 


2 


33 


2.5 


2.5 


2.5 


2 


3 


1.26 


0.112 


6 





2 


4 


0.124 


152 Gd 


2 


19 


2 


13 


1.8 


1.8 


3.2 


2 


1 


0.00 


0.104 


16 


10 


7 


19 


0.635 


154 Dy 


2 


23 


2 


15 


2.0 


2.0 


3.1 


2 


1 


0.75 


0.083 


26 


10 


7 


24 


0.530 


156 Er 


2 


32 


2 


25 


2.7 


2.8 


2.7 


2 


3 


1.24 


0.043 


20 


4 


5 


16 


0.450 



17 



TABLE I: (continued) 



nucleus R 4 / 2 


Ri/2 


Ro/2 


Ro/2 


R2/2 


-R2/2 /So a 


L 9 


Lp L 1 


n a 




exp 


th 


exp 


th 


exp 


th 










18bp t 


2.56 


2.42 


2.5 


3.7 


3.2 


2.4 1.71 0.085 


26 


6 


10 


25 0.813 


188p t 


2.53 


2.37 


3.0 


3.3 


2.3 


2.4 1.52 0.076 


16 


2 


4 


12 0.637 


190p t 


2.49 


2.28 


3.1 


3.4 


2.0 


2.3 1.42 0.015 


18 


2 


6 


15 0.637 


192p t 


2.48 


2.34 


3.8 


3.7 


1.9 


2.3 1.56 0.032 


10 





8 


12 0.681 


194p t 


2.47 


2.36 


3.9 


3.6 


1.9 


2.4 1.55 0.049 


10 


4 


5 


11 0.667 


196p t 


2.47 


2.33 


3.2 


2.9 


1.9 


2.3 1.37 0.079 


10 


2 


6 


11 0.639 


198p t 


2.42 


2.21 


2.2 


2.2 


1.9 


2.2 0.96 0.089 


6 


2 


4 


7 0.370 


200p t 


2.35 


2.00 


2.4 


2.0 


1.8 


2.0 0.00 0.000 


4 





4 


5 0.392 
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TABLE II: Same as Table 1, but for axially symmetric prolate deformed rare earth and actinide nuclei with R4/2 > 2.9 . Po, 
a, an d c are free parameters, related t o th e Davidson potential [Eq. (58 1], to the dependence of the mass on the deformation 
[Eq. ( 59 1] , and to the 7-potential [Eq. ( 40 1] . The theoretical predictions are obtained from the equations mentioned in subsec. 
|XI B| where further discussion can be found. 



nucleus 7? 


1/2 


R 


1/2 


Ro/2 


Ro/2 


R2/2 


R2/2 


Po 


c 


a 


L 3 


Lp 




n 


a 




exp 




th 


exp 


th 


exp 


th 


















15U Nd 


2.93 


3.13 


5.2 


7.9 


8.2 


5.8 


0.0 


2.1 


0.003 


14 


6 


4 


13 


2.012 


152 Sm 


3.01 


3.14 


5.6 


8.4 


8.9 


6.5 


0.0 


2.4 


0.000 


16 


14 


9 


23 


3.327 


154 Sm 


3.25 


3.27 


13.4 


13.0 


17.6 


18.6 


1.30 


6.9 


0.021 


16 


6 


7 


17 


0.515 


154 Gd 


3.02 


3.09 


5.5 


6.5 


8.1 


4.1 


0.0 


1.4 


0.024 


26 


26 


7 


32 


3.546 


156 Gd 


3 


24 


3 


25 


11.8 


10.8 


13.0 


14.3 


0.0 


5.3 


0.026 


26 


12 


16 


34 


0.933 


158 Gd 


3 


29 


3 


29 


15.0 


14.5 


14.9 


15.1 


1.99 


5.3 


0.025 


12 


6 


6 


14 


0.323 


160 Gd 


3 


30 


3 


30 


17.6 


17.3 


13.1 


13.2 


2.38 


4.5 


0.020 


16 


4 


8 


17 


0.125 


162 Gd 


3 


29 


3 


30 


19.8 


19.8 


12.0 


12.1 


2.52 


4.1 


0.008 


14 





4 


10 


0.078 


156 Dy 


2 


93 


3 


13 


4.9 


7.4 


6.5 


5.3 


0.0 


1.9 


0.014 


28 


10 


13 


31 


1.789 


158 Dy 


3 


21 


3 


22 


10.0 


9.6 


9.6 


10.3 


0.26 


3.8 


0.023 


28 


8 


8 


25 


0.496 


160 Dy 


3 


27 


3 


27 


14.7 


14.7 


11.1 


12.1 


1.92 


4.3 


0.005 


28 


4 


23 


38 


0.510 


162 Dy 


3 


29 


3 


30 


17.3 


15.7 


11.0 


11.2 


2.23 


3.8 


0.020 


18 


8 


14 


26 


0.742 


164 Dy 


3 


30 


3 


30 


22.6 


22.5 


10.4 


10.2 


2.68 


3.4 


0.000 


20 





10 


19 


0.100 


166 Dy 


3 


31 


3 


31 


15.0 


14.9 


11.2 


11.2 


2.39 


3.7 


0.047 


6 


2 


5 


8 


0.077 


160 Er 


3 


10 


3 


16 


7.1 


8.1 


6.8 


6.6 


0.00 


2.4 


0.013 


26 


2 


5 


18 


0.699 


162 Er 


3 


23 


3 


23 


10.7 


10.7 


8.8 


10.1 


1.29 


3.7 


0.013 


20 


4 


12 


23 


0.770 


164 Er 


3 


28 


3 


27 


13.6 


12.2 


9.4 


9.6 


1.83 


3.3 


0.026 


22 


10 


18 


33 


0.918 


166 Er 


3 


29 


3 


28 


18.1 


16.8 


9.8 


9.9 


2.22 


3.4 


0.002 


16 


10 


14 


26 


0.698 


168 Er 


3 


31 


3 


31 


15.3 


14.4 


10.3 


10.2 


2.29 


3.4 


0.041 


18 


6 


8 


19 


0.404 


170 Er 


3 


31 


3 


30 


11.3 


10.1 


11.9 


12.9 


1.64 


4.4 


0.083 


24 


10 


19 


35 


0.837 


i62 Yb 


2 


92 


3.07 


3.6 


6.8 


4.8 


4.0 


0.00 


1.4 


0.003 


24 





4 


15 


1.036 


i64 Yb 


3 


13 


3.18 


7.9 


8.3 


7.0 


7.4 


0.00 


2.7 


0.023 


18 





5 


13 


0.357 


166y b 


3 


23 


3.23 


10.2 


8.9 


9.1 


9.7 


0.66 


3.5 


0.038 


24 


10 


13 


29 


0.973 


168y b 


3 


27 


3.26 


13.2 


11.2 


11.2 


11.5 


1.52 


4.1 


0.028 


34 


4 


7 


25 


1.070 


170 Yb 


3 


29 


3.27 


12.7 


11.2 


13.6 


14.1 


1.36 


5.1 


0.035 


20 


10 


17 


31 


0.963 


172 Yb 


3 


31 


3.30 


13.2 


12.2 


18.6 


18.9 


1.66 


6.6 


0.055 


16 


10 


5 


17 


0.742 


174 Yb 


3 


31 


3.31 


19.4 


19.3 


21.4 


21.5 


2.44 


7.5 


0.019 


20 


4 


5 


16 


0.104 


176 Yb 


3 


31 


3.30 


13.9 


13.7 


15.4 


15.5 


1.97 


5.4 


0.036 


20 


2 


5 


15 


0.287 


178 Yb 


3 


31 


3.27 


15.7 


15.5 


14.5 


14.6 


1.88 


5.3 


0.000 


6 


4 


2 


6 


0.127 


166 Hf 


2 


97 


3.08 


4.4 


6.9 


5.1 


4.3 


0.00 


1.5 


0.006 


22 





3 


13 


0.873 


168 Hf 


3 


11 


3 


17 


7.6 


8.1 


7.1 


6.9 


0.00 


2.5 


0.023 


22 


4 


4 


16 


0.494 


170 Hf 


3 


19 


3 


21 


8.7 


8.7 


9.5 


8.8 


0.00 


3.2 


0.033 


34 


4 


4 


22 


0.970 


172 Hf 


3 


25 


3 


24 


9.2 


9.8 


11.3 


11.7 


0.00 


4.3 


0.031 


38 


4 


6 


26 


0.549 


174 Hf 


3 


27 


3 


25 


9.1 


10.4 


13.5 


13.6 


0.00 


5.0 


0.033 


26 


4 


5 


19 


0.832 


176 Hf 


3 


28 


3 


28 


13.0 


11.5 


15.2 


16.1 


1.31 


5.8 


0.038 


18 


10 


8 


21 


0.950 


178 Hf 


3 


29 


3 


28 


12.9 


12.3 


12.6 


13.0 


1.70 


4.6 


0.028 


18 


6 


6 


17 


0.356 


180 Hf 


3 


31 


3 


30 


11.8 


11.5 


12.9 


13.0 


1.92 


4.4 


0.068 


12 


4 


5 


12 


0.157 


176 w 


3 


22 


3 


21 


7.8 


9.1 


9.6 


9.5 


0.00 


3.5 


0.027 


22 


4 


5 


17 


0.881 


1,8 W 


3 


24 


3 


22 


9.4 


8.6 


10.5 


8.9 


0.00 


3.2 


0.039 


18 


10 


2 


15 


0.987 


1S0 W 


3 


26 


3 


25 


14.6 


13.1 


10.8 


11.5 


1.64 


4.2 


0.000 


24 





7 


18 


0.603 


182 W 


3 


29 


3 


29 


11.3 


11.5 


12.2 


12.5 


1.77 


4.3 


0.050 


18 


4 


6 


16 


0.195 


184 w 


3 


27 


3 


28 


9.0 


8.9 


8.1 


8.0 


1.57 


2.7 


0.080 


10 


4 


6 


12 


0.093 


186 W 


3 


23 


3 


25 


7.2 


7.2 


6.0 


6.3 


1.20 


2.1 


0.099 


14 


4 


6 


14 


0.130 


176 0s 


2 


93 


3 


10 


4.5 


6.9 


6.4 


4.6 


0.00 


1.6 


0.016 


24 


6 


5 


19 


1.747 


178 0s 


3.02 


3 


12 


4.9 


7.2 


6.6 


5.1 


0.00 


1.8 


0.017 


16 


6 


5 


15 


1.836 


180 Os 


3.09 


3 


22 


5.6 


7.1 


6.6 


6.9 


0.00 


2.4 


0.078 


10 


6 


7 


14 


1.021 


184 0s 


3.20 


3 


21 


8.7 


9.9 


7.9 


8.5 


1.21 


3.1 


0.011 


22 





6 


16 


0.886 


186 Os 


3.17 


3 


19 


7.7 


7.0 


5.6 


6.0 


0.00 


2.1 


0.063 


14 


10 


13 


24 


0.702 


188 Os 


3.08 


3 


15 


7.0 


7.2 


4.1 


4.4 


1.07 


1.5 


0.033 


12 


2 


7 


13 


0.170 


190 Os 


2.93 


3.07 


4.9 


5.6 


3.0 


3.1 


0.00 


1.0 


0.051 


10 


2 


6 


11 


0.419 
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TABLE II: (continued) 



nucleus R4./2 


-R4/2 


Ro/2 


Ro/2 


R2/2 


R2/2 


Po 


c 


L a 


Lp 




n a 




exp 


th 


exp 


th 


exp 


th 














228 Ra 


3.21 


3.24 


11.3 


11.0 


13.3 


13.3 


0.57 


5.0 0.016 


22 


4 


3 


15 0.177 




3.24 


3.26 


14.4 


14.3 


16.8 


17.0 


1.50 


6.4 0.002 


18 


2 


5 


14 0.214 


230r™ 


3.27 


3.27 


11.9 


11.6 


14.7 


14.7 


1.44 


5.3 0.034 


24 


4 


4 


17 0.243 


232rjij^ 


3.28 


3.28 


14.8 


14.0 


15.9 


16.5 


1.80 


5.9 0.022 


30 


10 


12 


31 0.426 


232-q 


3.29 


3.29 


14.5 


13.8 


18.2 


18.4 


1.74 


6.6 0.028 


20 


10 


4 


18 0.394 


234 u 


3.30 


3.30 


18.6 


18.3 


21.3 


21.8 


2.19 


7.8 0.011 


28 


8 


7 24 0.244 


236 jj 


3.30 


3.30 


20.3 


20.0 


21.2 


21.2 


2.38 


7.5 0.009 


30 


4 


5 


21 0.143 


238 jj 


3.30 


3.31 


20.6 


20.6 


23.6 


24.7 


2.38 


8.8 0.009 


30 


4 


27 43 0.665 


238p u 


3.31 


3.31 


21.4 


21.4 


23.3 


23.3 


2.61 


8.1 0.016 


26 


2 


4 


17 0.067 


240p u 


3.31 


3.31 


20.1 


19.9 


26.6 


26.6 


2.40 


9.4 0.018 


26 


4 


4 


18 0.117 


242p u 


3.31 


3.31 


21.5 


21.4 


24.7 


24.7 


2.52 


8.7 0.012 


26 


2 


2 


15 0.107 


248 Cm 


3.31 


3.31 


25.0 


24.8 


24.2 


24.3 


2.72 


8.5 0.004 


28 


4 


2 


17 0.159 


250 Cf 


3.32 


3.31 


27.0 


26.9 


24.2 


24.2 


2.88 


8.4 0.003 


8 


2 


4 


8 0.053 



TABLE III: Normalized [to the energy of the first excited 
state, E(2f)] energy levels of the ground state band (gsb) 
and the /3i and 71 bands of 162 Dy and 238 U, obtained from 
the Bohr Hamiltonian with /3-dependent mass for axially sym- 
metric prolate deformed nuclei using the parameters given in 
Table 2, compared to experimental data [3D]. See subsec. 
IXI Bl for further discussion. 





^Dy 


ib2 Dy 


238TJ 


238yj 




^Dy 


iW Dy 


238-q 


238yj 


L 


exp 


th 


exp 


th 


L 


exp 


th 


exp 


th 




gsb 


gsb 


gsb 


gsb 




7i 


7i 


7i 


7i 





0.00 


0.00 


0.00 


0.00 


2 


11.0 


11.2 


23.6 


24.7 


2 


1.00 


1.00 


1.00 


1.00 


3 


11.9 


12.1 


24.6 


25.5 


4 


3.29 


3.30 


3.30 


3.31 


4 


13.2 


13.3 


25.9 


26.7 


6 


6.80 


6.80 


6.84 


6.86 


5 


14.7 


14.7 


27.4 


28.1 


8 


11.41 


11.41 


11.54 


11.57 


6 


16.4 


16.5 


29.2 


29.8 


10 


17.04 


17.01 


17.27 


17.33 


7 


18.5 


18.5 


31.2 


31.7 


12 


23.57 


23.49 


23.97 


24.06 


8 


20.7 


20.8 


33.5 


33.9 


14 


30.90 


30.74 


31.51 


31.63 


9 


23.3 


23.3 


36.0 


36.3 


16 


38.90 


38.70 


39.82 


39.97 


10 


25.9 


26.0 


38.8 


39.0 


18 


47.58 


47.28 


48.78 


48.98 


11 


29.0 


28.9 


41.7 


41.9 


20 






58.31 


58.61 


12 


31.4 


32.1 


44.9 


45.0 


22 






68.31 


68.77 


13 


35.5 


35.5 


48.3 


48.3 


24 






78.71 


79.44 


14 


39.4 


39.9 


51.9 


51.8 


26 






89.46 


90.55 


15 






55.7 


55.5 


28 






100.57 


102.08 


16 






59.7 


59.4 


30 






112.10 


113.99 


17 






63.9 


63.4 












18 






68.2 


67.7 




Pi 


Pi 


Pi 


Pi 


19 






72.7 


72.0 





17.3 


15.7 


20.6 


20.6 


20 






77.3 


76.6 


2 


18.0 


16.7 


21.5 


21.6 


21 






82.1 


81.3 


4 


19.5 


19.0 


23.5 


24.0 


22 






87.0 


86.1 


6 


21.9 


22.6 






23 






91.9 


91.0 


8 


24.6 


27.4 






24 






97.0 


96.1 












25 






102.1 


101.3 












26 






107.4 


106.6 












27 






112.7 


112.0 
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TABLE IV: Comparison of experimental data [40] (upper line) for several B(E2) ratios of 7- unstable nuclei to predictions 
(lower line) by the Bohr Hamiltonian with /3-dependent mass (with 8 = A = 0), for the parameter values shown in Table 1. See 
subsec. IXI CI for further discussion. 



nucl. 






W, — >.fji 


1U1 — ^81 




2 2 — 


O2— >2i 




2 1 ->0 1 




2l-»0i 




2l-M)i 


x 10 3 


2i-Hh 


2i-H), 
X 10 3 
















a8 Ru 


1.44(25) 








1.62(61) 


36.0(152) 








1.82 


2.62 


3.42 


4.22 


1.82 


0.0 


1.36 


3.60 


100 Ru 


1.45(13) 








0.64(12) 


41.1(52) 


0.98(15) 






1.72 


2.40 


3.07 


3.73 


1.72 


0.0 


1.05 


10.89 


102 Ru 


1.50(24) 








0.62(7) 


24.8(7) 


0.80(14) 






1.78 


2.54 


3.28 


4.01 


1.78 


0.0 


1.27 


8.70 


104 Ru 


1.18(28) 








0.63(15) 


35.0(84) 


0.42(7) 






1.63 


2.18 


2.71 


3.21 


1.63 


0.0 


0.79 


22.41 


io 2pd 


1.56(19) 








0.46(9) 


128.8(735) 








1.76 


2.49 


3.19 


3.87 


1.76 


0.0 


1.22 


12.34 


104p d 


1.36(27) 








0.61(8) 


33.3(74) 








1.74 


2.45 


3.15 


3.85 


1.74 


0.0 


1.11 


8.13 


106p d 


1.63(28) 








0.98(12) 


26.2(31) 


0.67(18) 






1.85 


2.67 


3.49 


4.28 


1.85 


0.0 


1.49 


5.98 


108 Pd 


1.47(20) 


2.16(28) 


2.99(48) 




1.43(14) 


16.6(18) 


1.05(13) 


1.90(29) 




1.75 


2.45 


3.12 


3.75 


1.75 


0.0 


1.20 


15.82 


no pd 


1.71(34) 








0.98(24) 


14.1(22) 


0.64(10) 






1.76 


2.43 


3.01 


3.51 


1.76 


0.0 


1.31 


26.24 


106 Cd 


1.78(25) 








0.43(12) 


93.0(127) 








1.68 


2.32 


2.95 


3.58 


1.68 


0.0 


0.92 


10.44 


108 Cd 


1.54(24) 








0.64(20) 


67.7(120) 








1.85 


2.69 


3.52 


4.35 


1.85 


0.0 


1.49 


4.06 


110 Cd 


1.68(24) 








1.09(19) 


48.9(78) 




9.85(595) 




1.99 


2.97 


3.93 


4.87 


1.99 


0.0 


1.98 


1.61 


112 Cd 


2.02(22) 








0.50(10) 


19.9(35) 


1.69(48) 


11.26(210) 




2.00 


2.99 


3.98 


4.96 


2.00 


0.0 


1.99 


0.48 


114 Cd 


1.99(25) 


3.83(72) 


2.73(97) 




0.71(24) 


15.4(29) 


0.88(11) 


10.61(193) 




2.00 


2.99 


3.97 


4.94 


2.00 


0.0 


1.99 


0.74 


116 Cd 


1.70(52) 








0.63(46) 


32.8(86) 


0.02 






1.74 


2.46 


3.17 


3.90 


1.74 


0.0 


1.11 


4.42 


118 Cd 


>1.85 












0.16(4) 






1.71 


2.39 


3.06 


3.74 


1.71 


0.0 


1.00 


5.88 


118 Xe 


1.11(7) 


0.88(27) 


0.49(20) 


>0.73 












1.67 


2.28 


2.85 


3.39 


1.67 


0.0 


0.95 


21.93 


120 Xe 


1.16(14) 


1.17(24) 


0.96(22) 


0.91(19) 












1.60 


2.11 


2.60 


3.08 


1.60 


0.0 


0.67 


21.51 


122 Xe 


1.47(38) 


0.89(26) 


>0.44 














1.58 


2.05 


2.48 


2.89 


1.58 


0.0 


0.63 


29.29 


124 Xe 


1.34(24) 


1.59(71) 


0.63(29) 


0.29(8) 


0.70(19) 


15.9(46) 








1.59 


2.09 


2.57 


3.04 


1.59 


0.0 


0.63 


20.14 


128 Xe 


1.47(20) 


1.94(26) 


2.39(40) 


2.74(114) 1.19(19) 


15.9(23) 








1.63 


2.20 


2.75 


3.31 


1.63 


0.0 


0.73 


9.64 


1 "39 -c r 

" 2 Xe 


1.24(18) 








1.77(29) 


3.4(7) 








2.00 


3.00 


4.00 


5.00 


2.00 


0.0 


2.00 


0.00 


130 Ba 


1.36(6) 


1.62(15) 


1.55(56) 


0.93(15) 












1.56 


2.01 


2.41 


2.77 


1.56 


0.0 


0.61 


34.54 


132 Ba 










3.35(64) 


90.7(177) 








1.68 


2.30 


2.90 


3.50 


1.68 


0.0 


0.92 


15.21 


134 Ba 


1.55(21) 








2.17(69) 


12.5(41) 








1.75 


2.48 


3.21 


3.94 


1.75 


0.0 


1.14 


4.08 


142 Ba 


1.40(17) 


0.56(14) 
















1.55 


2.00 


2.41 


2.82 


1.55 


0.0 


0.49 


18.60 
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TABLE IV: (continued) 



nucl. 


4l^l 
2i-s-0i 


Bi-Mi 
2l->0i 


»i->t>i 

2l->0i 


lUl-J-8l 
2i^0i 




2l->0i 


2l->0, 

x 10 3 


U2-S-2! 
2i^0i 


2i->0i 

x 10 3 


I4 «Nd 


1.61(13) 


1.76(19) 









.25(4) 


9.3(17) 0.54(6) 


32.82(816) 




1.63 


2.17 


2.68 


3.15 


1 


.63 


0.0 


0.81 


26.86 


152 Gd 1.84(29) 


2.74(81) 









.23(4) 


4 2(8) 


2.47(78) 






1.98 


2.92 


3.81 


4.65 


1 


.98 


0.0 


1.95 


4.51 


154 Dy 1.62(35) 


2.05(42) 


2.27(62) 


1.86(69) 














1.91 


2.79 


3.64 


4.46 


1 


.91 


0.0 


1.70 


5.41 


156 Ej . 


1.78(16) 


1.89(36) 


0.76(20) 


0.88(22) 














1.70 


2.35 


3.00 


3.64 


1 


.70 


0.0 


0.98 


11.50 


192p t 


1.56(12) 


1.23(55) 






1 


.91(16) 


9.5(9) 








1.59 


2.09 


2.57 


3.05 


1 


.59 


0.0 


0.61 


16.98 


194p t 


1.73(13) 


1.36(45) 


1.02(30) 


0.69(19) 


1 


.81(25) 


5.9(9) 


0.01 






1.59 


2.09 


2.57 


3.04 


1 


.59 


0.0 


0.63 


19.78 


196p t 


1.48(3) 


1.80(23) 


1.92(23) 








0.4 


0.07(4) 


0.06(6) 




1.64 


2.21 


2.75 


3.28 


1 


.64 


0.0 


0.82 


20.83 


198p t 


1.19(13) >1.78 






1 


.16(23) 


1.2(4) 


0.81(22) 


1.56(126) 




1.82 


2.60 


3.36 


4.08 


1 


.82 


0.0 


1.41 


10.09 
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TABLE V: Comparison of experimental data [40] (upper line) for several B(E2) ratios of axially symmetric prolate deformed 
nuclei to predictions (lower line) by the Bohr Hamiltonian with /3-dependent mass (with 5 = A = 0), for the parameter values 
shown in Table 2. See subsec. IXI Dl for further discussion. 



nucl. 


4i — >2i 


6i 


«, 

°1 ' U l 


10i — >8i 


p L 


2 a — ► 2 I 

P L 


2a — »4i 
P L 


'2 — >(Ji 


2 7 — ^2i 




2i-+0i 


2j.-i.0i 


2i->0i 


2l-t-0i 


2H-01 
X 10^ 


X 10 J 


2 1 ->0, 

x 10^ 


2^0, 

x 10 


2 1 - > 0i 

x 10 J 


2 1 ->0, 

x 10 J 












lb4 Sm 


1.40(5) 


1.67(7) 


1.83(11) 


1.81(11) 


5.4(13) 




25(6) 


18.4(34) 




3.9(7) 




1.47 


1.69 


1.87 


2.06 


26.7 


50.0 


150 


47.5 


69.6 


3.7 


156 Gd 


1.41(5) 


1.58(6) 


1.71(10) 


1.68(9) 


3.4(3) 


18(2) 


22(2) 


25.0(15) 


38.7(24) 


4.1(3) 




1.48 


1.73 


1.95 


2.18 


29.7 


59.1 


191 


62.5 


92.4 


4.9 


158 Gd 


1.46(5) 




1.67(16) 


1.72(16) 


1.6(2) 


0.4(1) 


7.0(8) 


17.2(20) 


30.3(45) 


1.4(2) 




1.46 


1.66 


1.82 


1.98 


25.7 


45.9 


127 


64.0 


93.0 


4.8 


158 Dy 


1.45(10) 


1.86(12) 


1.86(38) 


1.75(28) 


12(3) 


19(4) 


66(16) 


32.2(78) 


103.8(258) 


11.5(48) 




1.50 


1.78 


2.04 


2.31 


30.5 


65.4 


232 


88.5 


131.7 


7.1 


160 Dy 


1.46(7) 


1.23(7) 


1.70(16) 


1.69(9) 


3.4(4) 




8.5(10) 


23.2(21) 


43.8(42) 


3.1(3) 




1.46 


1.68 


1.85 


2.03 


22.9 


43.5 


133 


78.6 


114.5 


6.0 


162 Dy 


1.45(7) 


1.51(10) 


1.74(10) 


1.76(13) 








0.12(1) 


0.20 


0.02 




1.45 


1.65 


1.80 


1.95 


23.9 


42.4 


116 


89.8 


129.8 


6.7 


164 Dy 


1.30(7) 


1.56(7) 


1.48(9) 


1.69(9) 








19.1(22) 


38.3(39) 


4.6(5) 




1.44 


1.62 


1.75 


1.86 


16.9 


29.1 


77 


99.7 


143.4 


7.3 


162 Er 










8(7) 




170(90) 


32.5(28) 


77.0(56) 


9.4(69) 




1.49 


1.75 


1.99 


2.24 


27.8 


58.3 


202 


91.1 


134.8 


7.2 


164 Er 


1.18(13) 




1.57(9) 


1.64(11) 








23.9(35) 


52.3(72) 


7.8(12) 




1.47 


1.70 


1.89 


2.09 


28.3 


53.5 


162 


103.8 


151.2 


7.9 


166 Er 


1.45(12) 


1.62(22) 


1.71(25) 


1.73(23) 








25.7(31) 


45.3(54) 


3.1(4) 




1.46 


1.66 


1.81 


1.96 


20.7 


38.2 


111 


100.0 


144.8 


7.4 


168 Er 


1.54(7) 


2.13(16) 


1.69(11) 


1.46(11) 








23.2(15) 


41.1(31) 


3.0(3) 




1.45 


1.65 


1.79 


1.93 


27.7 


47.2 


120 


100.6 


145.1 


7.4 


170 Er 






1.78(15) 


1.54(11) 


1.4(1) 


0.2(2) 


6.8(12) 


17.7(9) 




1.4(4) 




1.47 


1.69 


1.86 


2.03 


39.2 


67.9 


177 


78.6 


114.2 


5.9 


166 Yb 


1.43(9) 


1.53(10) 


1.70(18) 


1.61(80) 
















1.50 


1.78 


2.05 


2.33 


33.7 


71.0 


245 


97.2 


144.5 


7.8 


168 Yb 










8.6(9) 






22.0(55) 


45.9(73) 


8.6 




1.48 


1.72 


1.93 


2.14 


29.6 


57.5 


180 


82.9 


121.6 


6.4 


170 Yb 






1.79(16) 


1.77(14) 


5.4(10) 






13.4(34) 


23.9(57) 


2.4(6) 




1.47 


1.71 


1.91 


2.12 


30.6 


58.2 


176 


66.2 


97.1 


5.1 


172 Yb 


1.42(10) 


1.51(14) 


1.89(19) 


1.77(11) 


1.1(1) 


3.7(6) 


12(1) 


6.3(6) 




0.6(1) 




1.46 


1.67 


1.83 


1.99 


32.2 


55.9 


147 


51.6 


75.0 


3.9 


174 Yb 


1.39(7) 


1.84(26) 


1.93(12) 


1.67(12) 










12.4(29) 






1.45 


1.63 


1.75 


1.86 


20.9 


35.1 


88 


45.0 


64.9 


3.3 


176 Yb 


1.49(15) 


1.63(14) 


1.65(28) 


1.76(18) 








9.8 








1.46 


1.66 


1.82 


1.97 


27.9 


49.0 


132 


63.1 


91.6 


4.7 


174 Hf 










14(4) 




9(3) 


31.6(161) 


48.7(124) 






1.48 


1.74 


1.96 


2.20 


31.4 


62.2 


200 


66.9 


98.8 


5.3 


176 Hf 










5.4(11) 




31(6) 


21.3(26) 








1.47 


1.70 


1.89 


2.09 


30.8 


57.3 


169 


57.9 


84.9 


4.5 


178 Hf 




1.38(9) 


1.49(6) 


1.62(7) 


0.4(2) 




2.4(9) 


24.5(39) 


27.7(28) 


1.6(2) 




1.47 


1.69 


1.88 


2.07 


28.4 


53.1 


158 


73.8 


107.8 


5.6 


180 Hf 


1.48(20) 


1.41(15) 


1.61(26) 


1.55(10) 








24.5(47) 


32.9(56) 






1.46 


1.66 


1.82 


1.98 


34.9 


59.5 


151 


78.4 


113.4 


5.8 




1.43(8) 


1.46(16) 


1.53(14) 


1.48(14) 


6.6(6) 


4.6(6) 


13(1) 


24.8(12) 


49.2(24) 


0.2 




1.47 


1.69 


1.87 


2.04 


32.5 


58.3 


162 


79.9 


116.2 


6.0 


1W W 


1.35(12) 


1.54(9) 


2.00(18) 


2.45(51) 


1.8(3) 




24(3) 


37.1(28) 


70.6(51) 


4.0(4) 




1.48 


1.73 


1.95 


2.16 


40.7 


75.2 


216 


128.3 


187.3 


9.8 


"»W 


1.30(9) 


1.69(12) 


1.60(12) 


1.36(36) 








41.7(92) 


91.0(201) 






1.51 


1.80 


2.07 


2.34 


46.2 


91.9 


289 


165.7 


244.5 


12.9 


186 0s 


1.45(7) 


1.99(7) 


1.89(11) 


2.06(44) 








109.4(71) 


254.6(150) 13.0(47) 




1.53 


1.87 


2.20 


2.55 


39.7 


90.2 


335 


164.9 


247.4 


13.4 


188 0s 


1.68(11) 


1.75(11) 


2.04(15) 


2.38(32) 








63.3(92) 


202.5(304) 43.0(74) 




1.54 


1.89 


2.25 


2.63 


33.9 


83.9 


344 


229.8 


345.2 


18.7 
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TABLE V: (continued) 



nucl. 


<H2l 
2l->0i 


6i-»4i 
2l->0i 


8i->6i 
2n-0i 


10i->8i 
2i->0i 


2j3-»Ul 
2i->0i 

x 10 3 


2,3— >2x 
2i->0i 

x 10 3 


20 -t-»l 
2H-0, 

x 10 3 


21-1-0, 

x 10 3 


2 T ^2x 
2i-l0, 

x 10 3 


2 T ^4i 
21-1-0, 

x 10 3 




1.36(8) 








5.7(26) 




20(11) 15.6(59) 


28.1(100) 1.8(11) 




1.47 


1.70 


1.90 


2.09 


30.0 


56.4 


168 


63.6 


93.2 


4.9 


232 Th 


1.44(15) 


1.65(14) 


1.73(12) 


1.82(15) 


14(6) 


2.6(13) 


17(8) 


14.6(28) 


36.4(56) 


0.7 




1.46 


1.67 


1.84 


2.01 


25.8 


47.1 


135 


57.0 


83.0 


4.3 


234 U 
















12.5(27) 


21.1(44) 


1.2(3) 




1.45 


1.64 


1.78 


1.90 


20.7 


36.1 


97 


42.7 


61.8 


3.2 


236 


1.42(11) 


1.55(11) 


1.59(17) 


1.46(17) 
















1.45 


1.63 


1.76 


1.87 


19.3 


33.2 


87 


44.7 


64.5 


3.3 


238y 






1.45(23) 


1.71(22) 


1.4(6) 


3.6(14) 


12(5) 


10.8(8) 


18.9(17) 


1.2(1) 




1.45 


1.63 


1.75 


1.86 


18.9 


32.3 


83 


37.7 


54.5 


2.8 


238p u 










14(4) 




11(4) 










1.44 


1.62 


1.73 


1.84 


19.1 


31.7 


78 


41.6 


59.9 


3.0 


250 Qf 
















6.8(17) 


10.9(25) 


0.6(1) 




1.44 


1.61 


1.72 


1.81 


15.0 


24.9 


61 


40.0 


57.5 


2.9 



